%spanningen bepalen ten gevolge van een busfout
%for k=1:41
% voltages_lijn_1_fout = zeros(,30);
clear all
tic
einde_interval = 14;%# bussen bekijken
einde_interval2 = einde_interval;
teller = 0;
alle_spanningen_busfout = zeros(200*einde_interval,14);
tijd = zeros(200,einde_interval);
for k=1:einde_interval;
    
    filename = sprintf('busfout%03d',k);
    initpsat;
    Settings.plottype = 2;
    Settings.freq = 60;
    Settings.absvalues= 'on';
    runpsat(filename,'pf');
    Settings.fixt = 1;
    Settings.tstep = 1;%0.05; %TRAAAAAAAAAAAAAG
    
    Settings.freq = 60;
    Settings.absvalues= 'on';
    Settings.plottype = 2;
    runpsat('td');
    
    %begint vanaf index 0... (STOM!!!)
   v = genvarname('voltages_fout_bus', who);
   eval([v ' = Varout.vars']);%(:,(size(Bus.con,1)+1):2*size(Bus.con,1))']);
   lengte = size(Varout.vars,1);
   
   %spanning per tijdsstap per bus
   alle_spanningen_busfout(1+200*teller:200*teller+lengte,:) = Varout.vars(:,1:14);
   
   t = genvarname('tijd_fout_bus', who);
   eval([t ' = Varout.t']);%(:,(size(Bus.con,1)+1):2*size(Bus.con,1))']);
   
   %tijd per stap per lijn
   tijd(1:lengte,k) = Varout.t;
   closepsat;
   teller = teller+1;
   
   save busfout tijd alle_spanningen_busfout einde_interval2;
end
% initpsat;
% Settings.plottype = 2;
% Settings.freq = 60;
% Settings.absvalues= 'on';
% runpsat('IEEE_14_bus_original','pf');
% Settings.fixt = 1;
% Settings.tstep = 1;%0.05;
% Settings.freq = 60;
%     Settings.absvalues= 'on';
% runpsat('td');
% normaal = Varout.vars;
% toc